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ABSTRACT 

We deduce the equations that describe how polarized radiation is Comptonized by 
a hot electron gas. Low frequencies are considered, and the equations are expanded 
to second order in electron velocities. Induced scattering terms are included and a 
Maxwellian velocity distribution for the electrons is assumed. The special case of 
an axisymmetric radiation field is also considered, and the corresponding radiative 
transfer equations are found. Our results correct errors and misprints in previosly 
published transfer equations. The extension to a moving electron gas is made, and the 
radiative transfer equations are deduced to second order in gas velocity. We use the 
equations to study polarization in the Sunyaev-Zeldovich effect. 

Key words: polarization - radiative transfer - scattering - galaxies: clusters: general 
- cosmic microwave background 



1 INTRODUCTION 

A study of the spectral evolution of a radiation field caused 
by Compton scatterin g on an electron gas was first pub- 
lished by Kompaneets ( 1951: ) . He obtained the kinetic equa- 
tion describing the evolution of an isotropic and homoge- 
nous radiation field in an infinite, homogenous, thermal and 
nearly nonrelativistic electron gas. The Kompaneets equa- 
tion was later extended to include a non - homogenous radia- 
tion field in a pl ane parallel atmosphere ( Babuel-Peyrissac & 
[Rouvillois 196i: ). More recently the extension to arbitrarily 
high orders in frequency and temperature h as been discussed 
(^tebbins 1997|; [Challinor fc Lasenby 1998|). 



Compton scatteri ng of polarized radiation was stud- 
ied by Chan drasekhar (L960) in the zero temperature limit. 



Stark (1981) attempted to extend this work to nonzero tem- 
peratures. He obtained radiative transfer equations for the 
I and Q Stokes parameters for low frequency radiation scat- 
tered on a hot, nearly nonrelativistic thermal electron gas, 
including induced scattering and assuming an axisymmet- 
ric radiation field. The equations were derived by expanding 
to first order in kT/mc? and hv/m(? , where hv is photon 
energy, kT is thermal energy and mc? is the electron rest 
energy. 

The fully relativistic equations for all four Stokes pa- 
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rameters were obtained by Nagirner and Poutanen (199J) 
and fur ther ex tended to include induced scattering by Na- 
girner (1994a). These equations are very complicated and 
hard to use in practice. Therefore, when working with low 
frequencies and electron velocities, a simplified form is pre- 
ferred. Nagirner did also reduce the fully relat ivisti c equa- 
tions in this limit, to the same order as Stark (1981). 

However, when the simplified equations of Nagirner 



( 1994a ) are integrated assuming an axisymmetric field, the 
r esult is very different from the equations obtained by Stark 
(1981). The two sets of equations are supposed to describe 
the same phenomenon, but several terms disagree. To make 
the situation even more confusing , the eq uations in the Rus- 
sian original paper by Nagirner (|l994b|) are different from 
those in the English translation (Nagirner 1994a). As we 
show in this paper, the reason for the disagreements is that 
some relativistic corrections to the rotation angles between 
different polarization bases were neglected by Stark. Some 
small integration errors in both papers also contribute to 
the differences. 

In this paper we present the correct radiative transfer 
equations for polarized radiation in an electron gas to first 
order in hv/mc? and kT/mc? [equation (p^]. The equa- 
tions are derived from the Boltzmann coUisional equation. 
The special case of an axisymmetric field is considered [equa- 
tions dl^ ) and (|l5|)]. We then extend the equations to include 
a mo ving electron gas, to second order in gas velocity [equa- 
tion (^]. 

The transfer equations for Compton scattering of low 
frequency radiation on a relatively hot electron gas is of 
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high importance in cosmology. Clusters of galaxies contain 
a hot electron gas, and the cosmic background radiation is 
changed because of Compton scattering on this g as. This is 



the so called Sunyaev-Zeldovich e ffect (SZ effect) (Zeldovich 
& Sunjaev 1969; Rephaeli 1995). This effect together with 
X-ray data of clusters of galaxies is used e.g. to 
the Hubble constant. 



measure 

The SZ-effect for a moving electron 



gas ( Sunyaev fc Zeldovich 1980 ; Rephaeli 1995; Bazonov & 
Sunyae r 1998^ can be used to measure peculiar velocities 



of clusters. Such measurements have so far only used the 
intensity of the cosmic background radiation (CBR). Up- 
coming satellite experiments like MAP and Planck will try 
to measure polarization in the CBR. Therefore the equations 
describing the change in the Stokes parameters of the CBR 
caused by Compton scattering in clusters [equations ( |l3| ) 
and in this paper] will be important. We investigate 
if there is any temperature or velocity dependent polariza- 
tion produced. We also briefly discuss scattering of radiation 
from extended radio sources in the cluster on the intracluster 
gas and possible polarization produced from this effect. 

Here we derive the radiative transfer equations, start- 
ing in Section ^ with the Boltzmann coUisional equations 
in the rest frame of an electron. Using Lorentz transforma- 
tions, we transform the equations to the rest frame of the 
electron gas, and then expand to second order in electron 
velocity and first order in hv/mc?, before averaging over 
electron velocities (Section^. The final result is presented 
in equation ( p^ . We use different bases for the Stokes vector 
in the electron rest frame and in the rest frame of the gas. 
In Section ^ we discuss the transformation between these 
two bases. In Section |^, the equations are integrated in the 
case of an axisymmetric field to compare with the equa- 
tions of Stark equations ( p^ and (^^]. Then, in Section 
^ equation (^) is extended to a moving electron gas. We 
derive the radiative transfer equations to second order in 
gas velocity [equation (^^]. Finally, applications of the for- 
malism deduced in the previous sections to the SZ effect is 
discussed in Section ^. In the conclusions, we discuss the 
fiaws in the previously published results. Throughout this 
paper, the relativistic quantum system of units, in which 
the electron mass, speed of light and Planck's constant are 
given by m = c = ?i = 1, is used. 



2 THE BOLTZMANN COLLISIONAL 
EQUATION FOR THE STOKES 
PARAMETERS 

In this section, we study the evolution of a beam of polarized 
radiation in the 'system frame', in which the mean momen- 
tum of the electrons in the gas is zero. In this frame we 
introduce a fixed coordinate system {x, y, z) with arbitrary 
directions of the axes. A photon beam in the direction (6^, (p) 
(see Fig. |l|) is described by the 'number' Stokes vector 




n(i/, 0, 4>, t, x) 



where v is the frequency, x is position in space and t is time. 
The number Stokes parameters n, nq , nu and nv are related 
to the usual Stokes parameters by 




Figure 1. The incoming and outgoing photons in the xyz- 
coordinate system. 




Figure 2. The incoming photon beam and the electron velocity 
in the coordinate system where the z-axis is along the outgoing 
photon direction. 
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First we study the photon beam in the rest frame of 
an electron moving with velocity v. All quantites taken in 
this frame are denoted by subscript 'e'. The number Stokes 
parameters n,nQ ,nu ,nv are relativistically invariant, but 
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frequencies and angles transform according to the usual 
Lorentz transformations, 



7J/e(l + « cos Oe), COS O : 



COS ae + V 
1 +V COS Qfe 



■^"^{1 — V COS a)2 



COS 6e 



Ote 



(1 - cose) 



72(1 — 11 cos a){l ~ V cos q') ' 
1 / d 



7(1 — V cos q) 
7=(l-«^)-^/^ 



(1) 



where (a, <I>) are the polar and azimuthal angles between 
V and the photon [see Fig. |^ and dQ, is an infinitesimal 
solid angle. The direction unit vector uj gives the direction 
of the photon beam and Q is the angle between two photons 
which have angles a and a' with the electron velocity v. 
The Rnltzmann rolli.sirinal eqiiation describes the change of 



he{uc, ( L-d^.t^.Ti^^ caused bv scattering on particles with 



velocity v (Acquista & Anderson 1974). In the rest frame 
the equation takes the form 



dn ^ . 

— — + OJe ■ Ven = 

Ote 



ifie — dhe )Fe{v , Pe), 



(2) 



where and h~ are the number of photons per phase space 
volume scattered into the beam {ve,6e,4'e) and out of it, re- 
spectively from each scattering. Assuming an isotropic elec- 
tron gas, the number of electrons with velocity v is given 
by the electron distribution function Fe{v'^,pe), where pe 
is the electron density in the rest frame. Next, we assume 
that the distribution function can be written as Fe{v^ ,pe 



Pefe{v^), which in the system frame is F{v'^, p) = py ^f{v^ 
Here f{v'^) is normalised as 47r f{v^)v'^dv — 1. Going 



to 



the system frame, using the Lorentz transformations [equa- 
tion (^)] above, we find that 

dn 



dt 



+ u> ■ Vn 



{dUe — dne){'i — V cosa)pf{v ). 



(3) 



The rate of scattering from one photon beam (u'e, Q'e^'t'e) 
into {ve,Oe, cf>e) is givcu by 

dn^uldue = (1 -f N)LiReL2!^e ^ du'eii'edVL'e, 

where h' = h(v'9'4)'), dQ.L — dcosOLdSL and Re is the 



scatteri ng matrix in the re st frame, given by ( Berestetskii 
LifshitzpT Pitaevskii 197l|) 



3 (UeV- 

= 16^^"^ We) 



( Ru 

Rl2 



V 



Rl2 








R22 











Rs3 











R44 



where 



Ve 



16^ 

■7. ^e, 



R22 — cos 9e -I- 1, 



R\2 = cos Oe — 

Riz = 2 cos Qe 



(4) 



and 
R44 = cos Qe 



+ 



The Klein-Nishina cross section Re is given by 
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Re 



16tt 



(Jt 



— + COS^ Qe 



(5) 



where ar is the Thompson cross section. The L matrices 
transform between different polarization bases as discussed 
in the next section. The more photons there are in the state 
{i>e,de,<l>e), the morc likely it is that photons are scattered 
into this state. This fact (induced scatte ring) is included by 
the factor (1 + IM), where N is given by (Nagirner 1994a) 




The incoming and outgoing frequencies are related by the 
usual formula for the Compton frequency shift. 



" 1 — I^e(l — COS Qe) 

The rate of scattering out of the beam {ue, de, (t>e) is (Na- 



girner 1994a) 



dhe 



{Rlefle + ViLlR'e\.2fle)dQ.'e, 



where n" = h{y" ,6' Note that the spontanous scat- 
tering term (the first term on the right hand side) in this 
equation contains the Klein-Nishina cross section and not 
the scatt ering matrix. This is becaus e of the isotropy of the 
medium (Acquista & Anderson 1974). The cross section R'e 
and the scattering matrix are given by equation ^ and 
(^), except that Ve and v'e are replaced by v" and Ve respec- 
tively. The frequencies v" and Ve are related by the Compton 
formula for the scattering Ve v", 
_ Ve 

1 -I- I^e(l — C0s6e) ' 

Inserting the expressions for dnZ and dnt into equation (^ , 
we find that 



dh 
'dt 



+ a; ■ Vn 



p(l — V cos a)f{v ) 



(1 



^^e ~l ~ Ml n' I 

— — n — iteTi — NLiReL2n 

dl/e 



N)LiReL2 



3 POLARIZATION BASES 

So far, the polarization vector basis has not been consid- 
ered. The scattering matrix Re describes the scattering of 
photons he in a certain polarization basis which we call the 
'scattering basis'. This basis is dependent on the photons 



and electrons involved in ea ch scattering event ( Berestetskii 



Lifshitz & Pitaevskii 1971), and is therefore unconvenient 
when considering several scatterings. We wish to describe 
the polarization relative to the fixed set of axes [x, y, z) in 
the system frame. When we use the vector h in this frame, 
we use a polarization basis which we call the 'system basis'. 
In this basis, the Stokes parameters are measured relative 
to the meridian planes. The matrix L2 rotates the vector h 
from the system basis to the scattering basis, so that the 
matrix Re can be applied. Then Li rotates back to the sys- 
tem basis ( Chandrasekhar 1960). A Stokes vector can be 
rotated from one basis to another with the rotation matrix 



( [chandrasekhar 196C ) 
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L(X) = 



10 

COS 2x sin 2% 

— sin 2x cos 2% 

1 



where x is the rotation angle of the polarization basis vec- 
tors. Using the notation of Chandrasekhar, we write 

LlReL2 = L(7r - i2e)ReL(-jle), 

where vr — i2e and — zie are the rotation angles in the rest 
frame of the electron. These angles must be transformed to 
the system frame and expanded to second order in electron 
velocity. We wish to express iie and i2e in terms of the ze- 
r oth o rder rotation angles ii and 12 used by Chandrasekhar 
( 1960 ). We separate the factors which are dependent on elec- 
tron directions (which we call x and y) from ii and 12 which 
are not. This enables us, as will be shown in the next sec- 
tion, to collect all dependency on electron directions (which 
we average over) in one matrix, making the integration in 
Section W eas ier. The Lorentz transform ed rotation angles 
are given by (Nagirner & Poutanen 1993): 



cos lie = a;cosii -|- j/smii, 

sin lie = a; sin ii — 3/ cos ii 

cosi2e = a; cos 12 -|- y sin 22, 

sini2e = a; sin 22 — y cos 12, 



y = ^7"Vl-C^ sin$, 

2,2 , 
X + y = 1. 

Here the following definitions were used: 



(7) 



272(1 -wC)(l-uC') 



= cos Q = COS a' , 
■q = cos 6 ri' — cos 9' , 

= m 



•-c^'), (8) 



where a' is the angle between the photon (0' and the 
electron velocity v. The angle between the two photons is 
and /i = cos 0. The angles ii and 12 ar e the zeroth order 
rotation angles given by Chandrasekhar (1960), 



cos ii — 



cos 12 



■q-ri 



77' - TJfi 



smzi 



smz2 



1 


— ry2 






1 









(9) 



4 INTEGRATION OVER ELECTRON 
VELOCITIES 

The next step is to average over alle possible electron di- 
rections by integrating the right side of equation (m) over ^ 



from to 1 and $ from to 2tt and divide by 27r. Before 
integrating, all terms are expanded to second order in v and 
first order in u. Expanding to second order in v is equivalent 
to first order in kT /nu? when using a Maxwellian distribu- 
tion function for the electrons. We expand fi and h" near 
u": 



V ~ v' and v : 

h(u',e',(i)') 

h{v" ,6' ,(t>') 



n + 

au 

~i , dn' „ 
au 



n 



2 du'^ 



2 du^ 



where we redefine fi = n{v, 9' , (j)'). Using the Lorentz trans- 
formation equations (^, v'^ ju^ and dv'^/du^ are transformed 
to give 



du'e 
dUe 

dn' 



(i-O 



7(1-0 



(i-O 7(i-0 + '^(i-m)' 



7^(1-0' 



(7(1-0-^^(1-^))^ 



72(1-0)^' 

In appendix ^ we have expanded all quantities we need in 
equation (^ to first order in 1/ (low frequencies) and second 
order in v (equivalent to first order in temperature). The 
quantities u = 2xy and w = {2x^ — 1) are expanded instead 
of X and y for reasons which will become obvious. 

The rotation matrices Li and L2 contain cosine and sine 
of the angles 2(i2e — tt) and 2(— iie). Again we wish to sep- 
arate the factors that are dependent on electron directions 
from those that are not. We use equations (Q) to write 

cos2(7r — i2e) = (2x^ — 1) cos 2(7r — 12) — 2a;2/ sin 2(7r — 12), 

sin2(7r — j2e) = (2x^ — 1) sin 2(7r — 12) + 2a;i/ cos 2(7r — 12), 

cos2(— iie) = (2a;^ — 1) cos 2(— ii) — 2a;y sin 2(— ii), 

sin 2(— iie ) = (2x^ - l)sin2(-ii) + 2a;ycos2(-ii). (10) 

Using these equations we can write LiReL2 as 

LiReLa = L(7r-i2)AL(-ji), 

where A is given by, 

IOtt 

/ -Rii -R12TO 

— (-R22 + i?33)wW 



-fJT 



R12W 

~Rl2U 









R12U 

{R22 + R33)UW 

Risw^ - R22V? 








Rii 



First of all, one can see that all terms proportional to sin $ 
can be omitted since they will disappear in the averaging 
over <1>. The L matrices are now independent of electron di- 
rections, so they can be treated as constants in the integra- 
tion. Next, we insert the expanded expressions into equation 
(^) and use the relation (which can be found from geometric 
considerations) 

C' = K+ VI-mVI-C^ cos*. (11) 

Finally, Alathematica^^^ is used to integrate over $ and 
Then we average over velocities, using a non relativistic 
Maxwellian distribution, 



/(«') = (27rfcr)-=*/'e-"'''^'='^, 



(12) 
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where k is the Bohzmann constant, and T is the temperature 
of the gas. Remembering that {v) = and (u^) = SfcT = 3r 
where 



{X) 



we find that 
dh 



dt 



+ • Vn = — (7tp(1 — 2v)h 



dv 



+ TV 



2, av 



-2rL(7r - ii2)CL(-ii)n' 

+2(1 - [i)vm.{'K - i2)BL(-ji) 



2n' + .(f^ 



(13) 



where the B and C matrices are defined in appendix ^ 

These are the radiative transfer equations for all four 
Stokes parameters in a rather hot thermal isotropic electron 
gas (first order in kT /ra(?) for low frequency (first order in 
hv/mc?) radiation in the rest frame of the electron gas. If the 
equation for the I-Stokes parameter is integrated assuming 
no dependency on angles in the radiation field h (isotropic 



radiati on field ) we recover the Kompaneets equation (Kom- 
paneetf 1956). 



5 THE RADIATIVE TRANSFER EQUATION 
FOR AN AXISYMMETRIC FIELD 

Now the field h is assumed homogenous and symmetric 
about the z-axis, so there is no (f> dependence. With this 
symmetry, there is no f/-polarization, so only the / and Q 
components of the Stokes vector are considered. We also 
change to the intensity parameters (/, Q) by means of the 
relations 

/ dn 1/9/3 



47r^3 ' dv 4:Tvu^ V dv v 

d^n 12 - 6 9/1 d^i 



dv^ 4:TTV^ Atvv^ dv Atvv^ dv^ 



Making the substitution z = cos (0 — (j)'), we define 



D 



dz 



^ J-i VT 



dz 



^ J-i VT^^ 

2 r dz 



[L(^-i2)CL(^ii)]„,, 



where [M]^; means the upper left minor of the matrix M. 
Now, using — rjrj' + y^l — rj^y^l — r/'^z to integrate over 
z, we find 



D 
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5-3772-^2(3-5772) (5772 -3) (V^ - 1) 

I (772-l)(57,'2-3) 5(l- 772)(l-,7'2) 



and 

Fii = l- 377'2 + 3772(377'2- 1) + 

277^77'(3 - 577'2) + 27777'(377'2 - 1), 

F12 = -2(77'2-l) (1-3772 -37777' + 577^')'), 
F21 = -2(772 - 1)(1- 37777' - 377'2 + 577'^ 77), 
F22 = (77' - 1)('?'^ - 1)(3 - 107777'). 



Finally we define (using the notation of Stark 1981) 
X„ = J Xrj'"drj', 
and 

where X = {I,Q) and 77 = (0,1,2,3). With this notation, 
the radiative transfer equations for / and Q take the form 

1^ = -I{l-2v) + ^{{3^v')Io + {l-3v'KQo^l2 

-Q2) + F[Gi] + 2r[(3772 - l)(/o - 277/1 - 3/2) 
+2(5772 - 3)77(/3 + Q3 - Qi) + 2(3772 - l)(Qo 

-«2)] + ^ [/(-I + ^§^)Gr + 0(-l + -§^)Gq]}, 

(14) 



dQ 

ds 



-Q(l -2v) + —{{l- rf){h + 3Qo - 3O2 - 3/2) 
16 

+F[Gq] + 2r[(l - 772)(-2/o + 6(77/1 + h) 
-3(Qo - Q2) + 1077((3i - h - Qz))] 



(15) 



where s — arpt, and 



Gi = {3 - v )Io + {1 ~ Sri'){Qo - I2 - Q2) 
Gq = 



-7?((5 - 3772)/i + (3 - 5772)(Qi - h 



3)), 



(l- 772)(/o + 3Qo-3g2-3/2 
-77(3/1 + 5Qi- 5(53 -5/3)). 



Equations (^^ a nd (F s j) are in Section ^ compared with the 
resuhs of Stark (|l98]F 



6 THE RADIATIVE TRANSFER EQUATIONS 
FOR COMPTON SCATTERING OF 
POLARIZED LOW FREQUENCY 
RADIATION ON A MOVING ELECTRON 
GAS 

In Sections ^ and ^ we transformed the Boltzmann coUisional 
equations from the rest frame of the electron [equation (^] 
to the rest frame of the electron gas. Now we obtain the 
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equations in a frame of reference where the electron gas is 
moving. In this frame the electrons have a common nonzero 
mean velocity. Therefore, we now have to make a transfor- 
mation from the rest frame of the electron to a frame where 
the electron has a velocity which is a sum of the thermal 
velocity (the velocity in the rest frame of the electron gas) 
and the velocity of the gas. 

We expand to second order in thermal velocity (v) and 
electron gas velocity {vg). To this order, relativistic addition 
of velocities reduces to vt = v + Vg where vt is the total 
electron velocity observed from the moving frame, v is the 
velocity in the rest frame of the gas (thermal velocity) and 
Vg is the velocity of the gas observed in the moving frame. 
Using spherical coordinates, we add each of the components 
to get 

Vt sin at cos <^t = v sin a cos $ + sin Qg cos $9 , 
Vt sin at sin "l>t = v sin a sin $ + sin ag sin $ ^ , 

ut cos at = II cos Q + cos ftg , (16) 

where vt — |vt|, Vg = |vg|. Here (at,$t) are the angles 
between the total velocity vector and the photon, (q, $) are 
the angles between the electron velocity and the photon in 
the rest frame of the gas and (ag, ^g) are the angles between 
the gas velocity and the photon. 

We follow the same steps that we used to obtain the 
equations in the rest frame of the electron gas. The differ- 
ence is that we now transform to a frame where the electron 
velocity is vt = v -I- Vg instad of v. S o, i n equation (Iffl and 



in the expanded quantities [equation ( Al)], we have to make 
the following replacements [using equations (p^j: 



v^y 1 — cos $ 

w -^/l^^^ sin <!> 



~ cos* + Ug \/l - Cl cos ^g 

- sin $ + Ug ^1 - sin $g 



where (g = cos ag . Making these replacements, we write the 
right hand side of equation (^) as 

Fo + Fl(Cg, <S>g)Vg + f 1 (C, + F2(Cg, $g)«g 

+F2(C,<E>)«' + nC,i>,Cs,*9Hs, 

where the F-functions are just collections of terms propor- 
tional to different orders of v. When averaging over electron 
velocities, the term F((", "1?, ^g, <l?g) disappears since (v) = 
by definition. This is important, because the terms that 
are left are just a sum of terms dependent on gas veloc- 
ity {vgjCg^^g) ^nd terms dependent on thermal velocity 
(u, ^, $). There are no cross terms dependent on both! So in 
the radiative transfer equations for a moving electron gas, 
the part dependent on temperature (thermal velocity) and 
the part dependent on gas velocity are independent. 

This simplifies the calculations since one can calculate 
the terms dependent on thermal velocity and the terms de- 
pendent on gas velocity separately. The terms containing 
the gas temperature can be found assuming that the gas 
velocity is zero, and similarly the terms dependent on gas 
velocity can be found assuming that the temperature is zero. 
Then the different terms can be added to get the complete 
radiative transfer equations for a moving electron gas with 
nonzero temperature. 

We now find the terms dependent on gas velocity assum- 
ing the temperature to be zero. Going back to the Boltz- 



mann collisional equation in the rest frame of an electron 
[equations (^], one has to transform all quantities to the 
frame of reference where the electron has a velocity Vg . This 
is exactly what was done previously, except that one has to 
replace (vX,^) with {vg, (g^^g) in all the expressions and 
expansions used, since the only velocity now is the gas ve- 
locity. 

Since there is no averaging over electron velocities and 
directions this time (the electron gas has a constant velocity 
which all the electrons are following), one can change to a 
coordinate system in which the expressions and calculations 
are more easy. Previously, a coordinate system (x, y, z) with 
arbitrary axes was introduced. The incoming and outgoing 
photon directions were described with polar and azimuthal 
angles (77 = cos 9, 0) and {rj' = cos9' , (ft') with respect to this 
coordinate system. But since the directions of the axes were 
arbitrary, one can chose the z-axis to be in the direction 
of the gas velocity. With this choice of z-axis, the angle 
77 = cos 6 between the z-axis and the photon is just the angle 
C,g — cos ag between the photon and electron. The azimuthal 
angles are still called (j) and 0' for the photons with frequency 
V, and ly' . With this choice of axes, the change 77 — * Cs ^^'^ 
rj' — > ^g is made everywhere in the equations. 

The correction factors w and u from the Lorentz trans- 
formations of the angles ii and 12 need special attention. 
These are expressed (see appendix ^) through the angle 
<1>. That was a convenient angle when there was an inte- 
gration over electron directions. Now however, there is only 
an integration over incoming photon directions (note that 
dQ' — dC^'gdcj}' with the current choice of axes), so the angle 
(0 — 0') is a better choice than <&. The following relations 
between "I> and (0 — (f)') can be found from geometrical con- 
siderations: 



cos <3? 



Cg - CsM 



sin $ = — 



A/2 

■^g 



-4>')- 



Since there is no integration over electron velocities, it 
is now easy to find the equations describing the radiative 
transfer of polarized radiation in a moving zero tempera- 
ture electron gas. All that is necessary is to put the ex- 
panded Lorentz transformed quantities (appendix ^ into 
the Boltzmann equation (^, and replace (77, ?;',$) in the 
way described above. The resulting equations are the ra- 
diative transfer equations for a moving zero temperature 
electron gas. Adding the temperature dependent terms from 
equation ( p^ ) (with a coordinate system where the z axis is 
chosen to be along the gas velocity direction) we find that 

— — h ■ Vn, = —(JTpil — v(^ — 2u)n + 
at 



arp 



J- y" dC #'|li([1 + 21.(1 -M)]B-2rC 

+ [G + (2C' - <:)B]v + [H + (2C' - C)G - (1 - 3C'' 
-F2CC')B]i^')L2n' + Li (^[{ly + 4r)(l - /i) + v{C' - C)]B 

+ (C' - C)[(3C' - C)B + G]vA L2^^ + - 0' 
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+2r(l - At) WU^U'^ + 2NLiBL2(l - [i)v{2n 



dv 



dn'. 



(17) 



where we have written ^ instead of (g for simphcity. Here 
the matrices G and H come from the expansion of A, A = 
B + Gv + Hv^. These matrices are given in appendix M. The 
L matrices are again Li = L(7r — ^2) and L2 = L(— ii), where 
the angles ii and 12 are given by equations (remember to 
let ?? ^ Cs s-nd rj' — > (g). This is the equations describing the 
change of a radiation field upon interaction with a hot (first 
order in kT/mc?), moving (second order in Vg/c) electron 
gas to first order in hv/mc?. To be able to perform the 
integral over incoming electron directions, one can use the 
relation = CgCg + ^/l - (I^Jl- Cameos {(p - 0')- 



7 APPLICATION TO THE SZ EFFECT 

One important application of the equations deduced here, is 
the interaction of the cosmic background radiation with the 
hot electron gas in clusters of galaxies. In this case, the fre- 
quency is so low that only terms to zeroth order in frequency 
need to be taken into account. For small changes in the ra- 
diation field, the initial field can be inserted for n on the 
right hand side in eq uation (^^ and integrated along a path 
through the cluster ( Zeldovich fc Sunyaev 1969 ). The initial 
radiation field is the cosmic background radiation which can 
be assumed to be isotropic and homogenous. In this case, 
the integration in equation ( p^ ) over incoming photons can 
easily be made using the techniques of Section H. 

First we assume that the gas velocity is zero, studying 
only temperature dependent effects. In this case the solution 
of equation (^^ for small changes in the radiation field gives 
the following change An to first order in optical depth: 



An = 



-2rC] L2no + 4r(l - /i)LiBL2i^ 



dnp 



-r(l - ^l)\.^BL2ly■ 



2 d'^fiQ 



where s — arp is the optical depth and J ds is an integral 
through the optical depth of the cluster. The initial radiation 
field no is given by 



no = 




This is a matrix equation equivalent to four scalar equations 
each describing the change in the Stokes parameters I,Q,U 
and V. Making the integrations, the first equat ion gives the 
thermal SZ effect (IZeldovich & Sunyaev 19691). The other 



equations all give zero polarization to first order in optical 
depth as was excpected because of the isotropy of the cluster 
gas and the radiation field. 

The equation can also be used to study the effect of 
the intracluster medium on radiation from extended radio 
sources within the cluster. The state of polarization of the 



radiation coming from radio sources will be changed upon 
scattering with electrons in the intracluster gas. The change 
in the degree of polarization will of course be proportional 
to the optical depth. It will also be dependent on the tem- 
perature of the cluster gas, as can be seen by conside ring 
the change of n from a single scattering [from equation ( p7[ ) 
considering only the inscattering term]: 



_3 
16 



dC' j d<^'|Li[B- 2rC]L2no 



+At{1 - ^l)UBL2V^ + r(l - ^l)\.^BL2v■ 



• d ho 



where no now is the initial spectrum of radiation coming 
from the radio source. This time no is not isotropic since the 
radiation is only coming from some angles. A more detailed 
study is necessary to determine if the change of polarization 
of radiation from a radio source can be used to determine 
density and temperature of the cluster gas. 

Finally we consider the kinematic SZ effect which is the 
change of the background radiation because of scattering on 
a moving cluster gas. This time we assume that the tem- 
perature of the gas is zero. Again, using equation (^^, we 
find to first order in optical depth using an initial Planck 
function: 



An 



Anq 



V xe 
c - 

- — -I 

10 C2 ' 



Cs)^, 



where x = v/t. The first equation is just the formula for the 
kinem atic SZ effect which was given by Sunyaev & Zeldovich 
without derivation. The second equation shows that 



(^98^ 



some polarization is produced from the kinematic effect. The 
degree of polarization is proportional to the tangential ve- 
locity component of the c luster . This result was also given 
by Sunyaev & Zeldovich (^98^, again without deductions. 
The second result can be used to determine the tangential 
velocity of clusters of galaxies, but at the moment the po- 
larization degree produced (about 10~*) is too small to be 
observed. 



8 CONCLUSION 

In equations (^) we have obtained the radiative transfer 
equations for the Stokes parameters to first order in hv/mc? 
and k T /mc? . Comp aring these to those obtained by Na- 
girner (1994a; 19941), some small differences are seen. Com- 
paring the radiative transfer equations for an axisymmet- 
ric fiel d [equ ations ( p^ ) and (p^] with those obtained by 
Stark (1981), we find several disagreements. However, set- 
ting Q = Q, studying only the intensity of radiation, the 
equations actually agree. This indicates that the difference 
has to do with polarization, and indeed, in the paper of 
Stark, we find that the zeroth order rotation angles ii and 12 
were used instead of iie and i2e. The relativistic corrections 
to these angles were neglected, which, as is seen from the 
Taylor expansion, is altering the integration substantially. 
Including these corrections, the equations in the paper of 
Stark become identical to those obtained here [equations 
( p^ and dli])], except for the induced scattering terms. By 
carefully performing the last integration in that paper, we 



© 1999 RAS, MNRAS 000, §-| 



8 Frode K. Hansen and Per B. Lilje 



find that the induced scattering terms also agree with what 
we have derived here. 

Ma nipulating the equations in the paper of Nagirner 



(1994a) before integration, it is possible to show that they 
agree with equation of this paper. So there has to be an 
error somewhere in his integration. And indeed, by correct- 
ing the misprints in the integration, we find full agreement 
with the results presented here. So our equations (nsl) cor- 



rect the errors of Stark (1981) and Nagirner (1994a; 1994b) 



We also transformed the equations to a frame where 
the electron gas is moving [equations (|l7|)]. In this frame 
the equations were derived to second order in gas velocity. 
The transfer equations were tested by showing that they give 
zero polarization produced from the thermal SZ effect to first 
order in optical depth. We also showed that polarization 
from the kinematic SZ effect is very small. The formalism 
developed here can probably be used to study polarization 
of radiation from radio sources in the cluster upon scattering 
with the intracluster medium. 



where u = 2xy and w = {2x^ — 1). 



APPENDIX B: MATRICES FOR THE 
RADIATIVE TRANSFER EQUATIONS 

The following matrices were used in equation ( [l^ and (|l7 



/ + 1 


- 1 










+ 1 














2^ 





V 








2^ 


/ Cll Cl2 










Cl2 C22 













C33 





) 


\ 





C44 





Cll ^-1 + 2^ + 3^*^ - 2^^ C12 = -2(m - ifiii + 1), 



C22 = 1 + 2^-^^' - 2^f, 



C33 = 2 2^ - 2At^ 



ACKNOWLEDGEMENTS 

We thank Mark Birkinshaw for suggesting the study of po- 
larization in the SZ-effect and Juri Poutanen for help and 
comments. We also thank the Research Council of Norway 
for a travel grant. 



APPENDIX A: TAYLOR EXPANSIONS 

The following quantities expanded to first order in v and 
second order in v are needed in Section ^: 
2 

^ l-2z.(l-/i). 



dl/e 



1 + 21.(1 -/i), 

Ki-m) + ^(C'-0 + ^\'(C'-C), 
-Ki-m) + ^'(C'-C) + ^'\'(C'-C), 



-—^ ^ ^2 ^ ^ (C - C) , 
§ - l + 2<-.^(l-3n, 
7 ~ 1 + -V , 

COS Be ^ li)(v{C + C,') 

W(C' + C" + CC'-i)), 

COS^ Be « ^JL^ - V2{1 - ^)^(C + C') 

-^'(1-m)(-2m-C' + 3K' 

-2CC' + 4KC'-C" + 3K"), 

« i + 2w^i— ^(C^ - l)(cos^$- 1) 
1 + M 



u ^ 2v{l ~ -^—^ sin $ + 



2v^^^J^^{C + C')sm$, (Al) 



H = 

Gil 
G13 

G23 

ii^ii 
H12 



Gil 


Gil 


G13 


Gil 


Gil 


'-'23 


— G13 


— G23 


G33 












v 

/ i?ll 

-ffl2 

^^13 —H23 His 



H12 His 
H22 H23 











V 

2{fl-l)^l{(:, + Q, 






G33 




H44, 



= 2ifi~i)^i^aV'i--c/ 



= 2(M + l)v/l-Clv/l-a' 



G33 



Gil 



(1 - + - m(-2 + 3C,' + 4C,C; + 3^')], 

+ (c. - c;)' + 2m(-i + ci - c.c; + c;') 
V(3c'+4c,c;+3c;')], 

H22 = -4: + 2^i + 2fi' + {5-A^i + 3^i^)((:g+C'g^) 
+ (2-14^* + 4/i')CsC;, 

2 



H33 = 

i?13 = 

H44 — 



^2-2^ + 4/i' + 2(1 + ^)iC', + ) - 2(1 + 3m)CsC;, 
1 + 2m, 



iCa+Cg)Gri, H23 ^ ^ ^ _ \ Cg + C'g)G23, 



2(i-m)(i-C9-C -CsC^). 



1 + fiV 1 — A* 
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